Alpha beta filter

An alpha beta filter (or alpha-beta filter) is a simplified form of observer for estimation, data smoothing and control applications. It is closely related to Kalman filters and to linear state observers used in control theory. Its principal advantage is that it does not require a detailed system model.

Contents

Filter equations

An alpha beta filter presumes that a system is adequately approximated by a model having two internal states, where the first state is obtained by integrating the value of the second state over time. Measured system output values correspond to observations of the first model state, plus disturbances. This very low order approximation is adequate for many simple systems, for example, mechanical systems where position is obtained as the time integral of velocity. Based on a mechanical system analogy, the two states can be called position x and velocity v. Assuming that velocity remains approximately constant over the small time interval ΔT between measurements, the position state is projected forward to predict its value at the next sampling time using equation 1.

\text{(1)} \quad \hat{\textbf{x}}_{k} \leftarrow \hat{\textbf{x}}_{k-1} %2B \Delta \textrm{T}\ \textbf{ } \hat{\textbf{v}}_{k-1}

Since velocity variable v is presumed constant, so its projected value at the next sampling time equals the current value.

\text{(2)} \quad \hat{\textbf{v}}_{k} \leftarrow \hat{\textbf{v}}_{k-1}

If additional information is known about how a driving function will change the v state during each time interval, equation 2 can be modified to include it.

The output measurement is expected to deviate from the prediction because of noise and dynamic effects not included in the simplified dynamic model. This prediction error r is also called the residual or innovation, based on statistical or Kalman filtering interpretations

{\textbf{(3)}} \quad \hat{\textbf{r}}_{k} \leftarrow \textbf{x}_{k} - \hat{\textbf{x}}_{k}

Suppose that residual r is positive. This could result because the previous x estimate was low, the previous v was low, or some combination of the two. The alpha beta filter takes selected alpha and beta constants (from which the filter gets its name), uses alpha times the deviation r to correct the position estimate, and uses beta times the deviation r to correct the velocity estimate. An extra ΔT factor conventionally serves to normalize magnitudes of the multipliers.

\textbf{(4)} \quad \hat{\textbf{x}}_{k} \leftarrow \hat{\textbf{x}}_{k} %2B (\alpha)\ \textbf{r}_{k}
\textbf{(5)} \quad \hat{\textbf{v}}_{k} \leftarrow \hat{\textbf{v}}_{k} %2B ( \beta / [ \Delta \textrm{T} ] )\ \textbf{r}_{k}

The corrections can be considered small steps along an estimate of the gradient direction. As these adjustments accumulate, error in the state estimates is reduced. For convergence and stability, the values of the alpha and beta multipliers should be positive and small.

\quad 0 < \alpha < 1
\quad 0 < \beta < 1

Values of alpha and beta typically are adjusted experimentally. In general, larger alpha and beta gains tend to produce faster response for tracking transient changes, while smaller alpha and beta gains reduce the level of noise in the state estimates. If a good balance between accurate tracking and noise reduction is found, and the algorithm is effective, filtered estimates are more accurate than the direct measurements. This motivates calling the alpha-beta process a filter.

Algorithm Summary

Initialize.

Update. Repeat for each time step ΔT:

  Project state estimates x and v using equations 1 and 2
  Obtain a current measurement of the output value
  Compute the residual r using equation 3
  Correct the state estimates using equations 4 and 5
  Send updated x and optionally v as the filter outputs

Sample Program

Alpha Beta filter can be implemented in C[1] as follows:

#include<stdio.h>
#include<stdlib.h>

main()
{
	float dt = 0.5;
	float xk_1 = 0, vk_1 = 0, a = 0.85, b = 0.005; 

	float xk, vk, rk;
	float xm;

	while( 1 )
	{
		xm = rand() % 100;// input signal

		xk = xk_1 + ( vk_1 * dt );
		vk = vk_1;

		rk = xm - xk;

		xk += a * rk;
		vk += ( b * rk ) / dt;

		xk_1 = xk;
		vk_1 = vk;

		printf( "%f \t %f\n", xm, xk_1 );
		sleep( 1 );
	}
}

Result

The corresponding images depicts the outcome of above program in graphical format. In first image, the blue worm is input signal and red is output; in second image, blue worm is input signal and yellow is output; and, in third image, the blue worm is input and green is output. For first two images, the output signal is visibly smoother than the input signal. Output avoids extreme spokes generated in input; rather following a more subtle path to have a similar effect. Also, the output moves in an estimate of gradient direction of input.

Higher is the alpha parameter, higher is the effect of input x and lesser is the dampening. Low value of beta is effective in controlling sudden surges in velocity. Also, as alpha increases beyond unity, we observe the opposite effect in output. The output gets rougher and uneven.[2]

Relationship to general state observers

More general state observers, such as the Luenberger observer for linear control systems, use a rigorous system model. Linear observers use a gain matrix to determine state estimate corrections from multiple deviations between measured variables and predicted outputs that are linear combinations of state variables. In the case of alpha beta filters, this gain matrix reduces to two terms. There is no general theory for determining the best observer gain terms, and typically gains are adjusted experimentally for both.

The linear Luenberger observer equations reduce to the alpha beta filter by applying the following specializations and simplifications.

Relationship to Kalman Filters

A Kalman filter estimates the values of state variables and corrects them in a manner similar to an alpha beta filter or a state observer. However, a Kalman filter does this in a much more formal and rigorous manner. The principal differences between Kalman filters and alpha beta filters are the following.

A Kalman filter designed to track a moving object using a constant-velocity target dynamics (process) model (i.e., constant velocity between measurement updates) with process noise covariance and measurement covariance held constant will converge to the same structure as an alpha-beta filter. However, a Kalman filter's gain is computed recursively at each time step using the assumed process and measurement error statistics, whereas the alpha-beta's gain is computed ad hoc.

The alpha beta gamma extension

It is sometimes useful to extend the assumptions of the alpha beta filter one level. The second state variable v is presumed to be obtained from integrating a third acceleration state, analogous to the way that the first state is obtained by integrating the second. An equation for the a state is added to the equation system. A third multiplier, gamma, is selected for applying corrections to the new a state estimates. This yields the alpha beta gamma update equations.

\hat{\textbf{x}}_{k} \leftarrow \hat{\textbf{x}}_{k} %2B ( \alpha )\ \textbf{r}_{k}
\hat{\textbf{v}}_{k} \leftarrow \hat{\textbf{v}}_{k} %2B ( \beta / [ \Delta \textrm{T} ] )\ \textbf{r}_{k}
\hat{\textbf{a}}_{k} \leftarrow \hat{\textbf{a}}_{k} %2B ( \gamma / [2\ \Delta \textrm{T} ^ \textrm{2}] )\ \textbf{r}_{k}

Similar extensions to additional higher orders are possible, but most systems of higher order tend to have significant interactions among the multiple states, so approximating the system dynamics as a simple integrator chain is less likely to prove useful.

See also

References

[3] [4] [5] [6] [7] [8] [9]

References

  1. ^ Tremor Cancellation in Handheld Microsurgical Devices, TC83 by Gaurav Mittal, Deepansh Sehgal and Harsimran Jit Singh, Punjab Engineering College
  2. ^ Tremor Cancellation in Handheld Microsurgical Devices, TC83 by Gaurav Mittal, Deepansh Sehgal and Harsimran Jit Singh, Punjab Engineering College
  3. ^ "The Alpha-Beta Filter", Penoyer, Robert; C Users Journal, July 1993
  4. ^ "Engine Speed Monitoring: The Alpha-Beta Filter", Microstar Laboratories
  5. ^ "Reconciling steady-state Kalman and alpha-beta filter design", by Painter, J.H.; Kerstetter, D.; Jowers, S. IEEE Transactions on Aerospace and Electronic Systems, Volume 26, Issue 6, Nov 1990 Page(s):986 - 991
  6. ^ Fixed-lag alpha-beta filter for target trajectory smoothing Ogle, T.L.; Blair, W.D. Aerospace and Electronic Systems, IEEE Transactions on Volume 40, Issue 4, Oct. 2004 Page(s): 1417 - 1421
  7. ^ "Description of an alpha-beta Filter in Cartesian Coordinates", by Cantrell, Ben H., NAVAL RESEARCH LAB WASHINGTON DC, March 21, 1973
  8. ^ "Comparison of Four Filtering Options for a Radar Tracking Problem", by Lawton, John A. ; Jesionowski, Robert J. ; Zarchan, Paul. NAVAL SURFACE WARFARE CENTER DAHLGREN DIV VA, 1979.
  9. ^ "Mathematical Techniques in Multisensor Data Fusion", By David Lee Hall, Sonya A. H., Artech House, 2004, ISBN 1-58053-335-3, section 4.4.4

External links